hep-ph/9809483 
RUB-TPII-14/98 



Hard exclusive electroproduction of two pions and their 

resonances 

M.V. Polyakov 

Petersburg Nuclear Physics Institute, Gatchina, St. Petersburg 188350, Russia 

and 

Institut fur Theoretische Physik II, Ruhr-Universitat Bochum, D-44780 Bochum, 

Germany 



Abstract 

We study the hard exclusive production of two pions in the virtual photon 
fragmentation region with various invariant masses including the resonance region. 
The amplitude is expressed in terms of two-pion light cone distribution amplitudes 
(27tDA's). We derive dispersion relations for these amplitudes, which enables us to 
fix them completely in terms of tttt scattering phases and a few low-energy subtrac- 
tion constants determind by the effective chiral lagrangian. Quantitative estimates 
of the resonance as well tttt background DA's at low normalization point are made. 
We also prove certain new soft pion theorem relating two-pion DA's to the one- 
pion DA. Crossing relations between 2-7rDA's and parton distributions in a pion are 
discussed. 

We demontrate that by studying the shape of the tttt mass spectra (not absolute 
cross section!) in a diffractive electroproduction one can extract the deviation of the 
meson (tt, p, etc.) wave functions from their asymptotic form 6z(l — z) and hence to 
get important information about the structure of mesons. Also we discuss how the 
moments of quark distributions in a pion can be measured in the hard diffraction. 
We suggest (alternative to Soding's) parametrization of tttt spectra which is suitable 
for large photon virtuality. 
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1 Introduction 



Recently it became possible to measure with high precision two pion hard diffractive pro- 
duction in the deeply virtual photon fragmentation region. The paper contains a detailed 
study of the "upper blob" of this process, more precisely, of the two-pion distribution 
amplitude. 

It has been realized recently [|1], |2|, |3| that a number of hard exclusive processes of the 
type: 

7l(<?) + Tip) — > F{q') + T'(p') , (1.1) 

at large Q 2 with t = [p — p') 2 , %Bj = Q 2 /2(pq0 and q' 2 = M F , fixed are amenable to QCD 
description. The factorization theorem |IJ asserts that the amplitude of the process (|1 . 1| ) 
has the following form: 

J dz J dxfJj' T (x,x - x B f,t,ii) H ij (Q 2 x/x B j,Q 2 ,z,fi) ®f (z,pi) 

+power-suppressed corrections, (1.2) 

where ff, T is T — > X" skewed parton distribution [|, |5|, [TJ, 0, ^ (for a review see ||), &f(z, fi) 
is the distribution amplitude of the hadronic state F (not necessarily one particle state), 
and Hij is a hard part computable in pQCD as series in a s (Q 2 ). Here we shall study 
general properties of the distribution amplitudes $J(z,[i) when the final hadronic state 
F is a two pion state (F = tttt). 

The 27rDA's were introduced recently in |7]] in the context of the QCD description 
of the process 7*7 — > 2tt, first model estimate of the corresponding DA at the two-pion 
invariant mass below the resonance region was done in . In the present paper we discuss 
how to describe the resonance region. We shall see that calculations of the two-pion DA in 
low energy region H and the knowledge of irn scattering phases are sufficient to determine 
27rDA's in a wide range of two-pion invariant masses. 

For definiteness we shall consider hard electroproduction of two pions off the nucleon: 

72 + A^2vr + T', (1.3) 

where the final state T' can be N, A, ttN, etc. The usage of 27rDA's to describe the process 
( |1.3|) has several advantages comparing to the case of single resonance DA's (p , p' , etc.): 

• The 27rDA's are the most general objects to describe the form of the tttt mass spectra 
in diffractive production of pions. We shall see that by studying the shape of the 
two-pion mass spectra (not the absolute cross section!) in diffractive production of 
two pions one can extract the information about the deviation of the meson (tt, p) 
DA's from their asymptotic form, i.e. essentially non-perturbative information about 
the structure of mesons. 
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• The usage of 27rDA's is anyhow necessary in studies of channels with strong nonres- 
onant background, for example, in the case of hard electroproduction of two neutral 
pions. 

• By crossing symmetry the 27rDA's are related to skewed parton distribution (SPD's) 
in a pion (see section Were 27rDA's measured one could get information about 
SPD's in a pion. 

• From theoretical point of view 27rDA's are necessary to define DA's of unstable 
resonances (p, p ; , f 2 , P3, etc.). 

The plan of this paper is as follows. Section ^| contains definitions and general prop- 
erties of two-pion distribution amplitudes. In Section |3| we prove soft pion theorems for 
the 27rDA's. These theorems relate the 27rDA's to the DA's of a single pion in a limit 
when one of the produced pions is soft. In section [| we derive relations between 27rDA's 
and quark distributions (skewed and not) in a pion. These relations can be, in principle, 
used for measurements of quark distributions in a pion in hard exclusive reactions. Also 
they constrain the 27rDA's by already known parton distributions in a pion. Section [| 
is devoted to derivation of the dispersion relations for the 27rDA's. These dispersion re- 
lations allow us to express the 27rDA's in a wide range of invariant two-pion masses in 
terms of (known) pion-pion phase shifts and a few low-energy subtraction constants. In 
subsection |5.1| we explain how to extract resonance DA's from 27rDA's. In section |6] we 
estimate the low-energy subtraction constants in the instanton model of the QCD vac- 
uum, also results for p-meson DA's are presented. In section [7| we discuss the applications 
and implications of the 27rDA's to the hard electroproduction of pions. Conclusions and 
an outlook are given in section []. 

2 Two-pion light-cone distribution amplitudes. 

We shall consider the following twist-2 chirally even and odd two-pion distribution am- 
plitudes (2ttDA's) 

<f>f(z,(,w 2 ) = —fdx-e-^ zP+x ~ 
11 An J 

x out (ir a ( Pl )iT b (p 2 )\^(x)nT^(0)\0}l +=x =q , (2.1) 



<n 6 (2,cv) 



Anw 2 



dx 



zP+3 



x out(n a ( Pl ) n\ P2 ) | 4>{x) vnTT^O) |0> 



=x±=0 



Here, n is a light-like vector 



n 



(2.2) 



0), which we choose to be n u = (1,0,0, 1). For 



3. 



any vector, V, the "plus" light-cone coordinate is defined as V + = {n ■ V) = V° + V 
the "minus" component as V~ = V° — V 3 . The outgoing pions have momenta Pi,p 2 , arid 



P = pi + p 2 is the total momentum of the final state, with P 



w 



Finally, T is a 
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flavor matrix (T — 1/2 for the isosinglet, T = r 3 /2 for the isovector 27rDA). [In the above 
equations the path-ordered exponential of the gauge field, required by gauge invariance 
is always assumed]. The constant of dimension of mass is defined as matrix element 
of local operator: 

lim (7r a ( Pl ) tt V) I HO) a* ^(0) |0) = e ab3 ^(p^ 2 - tftf). (2.3) 



The generalized distribution amplitudes, eqs. (|2.1|j2.2j ), depend on the following kine- 



matical variables: the quark momentum fraction with respect to the total momentum of 
the two-pion state, z; the variable £ = pf / P + characterizing the distribution of longitu- 
dinal momentum between the two pions, and the invariant mass (cm. energy), w 2 = P 2 . 
In what follows we shall work in the reference frame where Pi =0. In this frame 



'I' 2 _ W 2 (l-() 2 _ 2 2 



P~ = p+, Pi = p+ ~ , Pi = Pt± = P2± = ^C(l-C)-< (2-4) 
From these relations one obtains the following kinematical constraint: 

C(l-C)>^|- (2.5) 



The isospin decomposition of the 27rDA's eqs. Q2.1| , |2.2|) has the formQ: 



$ ab = 5 ab tr(T)$ /=0 + ^tr([r a , r 6 ]T)$ /=1 . (2.6) 
From the C-parity one can easily derive the following symmetry properties of the isoscalar 



(I = 0) and isovector (I = 1) parts of 2ttDA's eqs. (0,0) 



&= (zX,w 2 ) = ^ I = (l-z,C,w 2 ) = ^ I =°(z,l~C,w 2 ), 
&=\zX,w 2 ) = $ /=1 (l-z,CV) = -$ /=1 M-CV). (2.7) 



The isospin one parts of 27rDA's eqs. fl2.1| , |2.2|) are normalized as follows: 



£dztff\z,C,w 2 ) = (2{-l)F v (w 2 ), 

f 1 dz^r\z, C, w 2 ) = (2C - l)F t {w 2 ) , (2.8) 
Jo 

where F^iyo 2 ) is the pion e.m. form factor in the time-like region (^(0) = 1) and F t (w 2 ) 
is the "tensor" form factor of the pion normalized by F t (0) = 1. 

Let us decompose 27rDA's in eigenfunctions of the ERBL Q evolution equation (Gegen- 
bauer polynomials C^ 2 (2z — 1)) and in partial waves of produced pionsQ (Gegenbauer 

1 We do not specify the indices (|| and _L) if an equation holds for both chirally even and odd DA's 
2 In the formulae below we, for simplicity, neglect threshold effects. How to take them into account is 
discussed in the Appendix. 
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1/2 

polynomials C l ' (2( — 1), see below). Generically the decomposition (for both isoscalar 
and isovector DA's) has the form: 

oo n+1 

«&(*, C, w 2 ) = 6z(l - z) £ £ iU^Cf 2 ^ - 1)C, 1/2 (2C - 1). (2.9) 

n=0 Z=0 

Using the symmetry properties (|2.7f) we see that the index n goes over even (odd) and 
index I goes over odd (even) numbers for isovector (isoscalar) 27rDA's. The normalization 
conditions ( |2.8|) correspond to Bq^w 2 ) = F n (w 2 ) and -B^(O) = 1. 

The Gegenbauer polynomials C^J 2 (2z — 1) are eigenfunctions of the ERBL || evolution 
equation and hence the coefficients B n \ are renormalized multiplicatively (for even n and 
odd/): 



n , 2 x „ / 2 n V 7 ' l - 7o)/{2/3o) 



where /3o = 11 — 2/3 iVy and the one loop anomalous dimensions are JITJ 

2 



7n 



7n" 



3 v 



[l + n)(2 + n) 

n+1 i 



n+1 i 
fc=2 K 



k=2 



(2.10) 



From the decomposition fl2.9|) and eq. ( |2.10|) we can immediately make a simple pre- 
diction for the ratio of the hard P— and F— waves production amplitudes of pions in the 
reaction — > 27rp at large virtuality of the incident photon and fixed w 2 and Bjorken 



x: 



F— wave amplitude 1 

P-wave amplitude ~ l O g(g 2 ) 50 /("- 6iV /) ' ^' ^ 



or generically for the 2/c + 1 wave: 

(2k + 1)— wave amplitude 1 

P— wave amplitude log(Q 2 ) 7 2'=^ 2/3o ' ) 

3 Production of soft pions 



(2.12) 



In this section we prove the low-energy theorem for the isovector $y 1 (z,(,w 2 ) relating 
this 27rDA to the single pion DA (p w (z) [IT]. The latter is defined as: 



r 3 ,1 



(7r°(P)|^(x)n 7 5^(0)|0> = iUP + / dz e izP ^(z). (3.1) 
2 Jo 



Here f w = 93 MeV is the pion decay constant. 

Let us consider the matrix element out {^ a {pi) ^ b {p2) \^{.x) if)(0) |0) entering the 
definition of the chirally even isovector 27rDA eq. ( |2.1[ ) and consider the kinematical limit 
when the 4-momentum of one of the pions goes to zero. This limit corresponds to w 2 — > 
and C — ► 1 (or £ — > 0). Using the LSZ reduction formula we can write (for definiteness 
we take p% — » which corresponds to «j 2 -> and ( — > 1): 

3 

o^ a (pi) 7r 6 (p 2 ) | ^(x) n T - ^(0) |0) = z I d A ye i{ v^\dl + m*) 

x (n^p^T^M^x) (3-2) 

where <p % {y) is a pion interpolating field which, owing to PCAC, can be related to the 
divergence of the axial current <j)\{y) = f ^ 2 9 ll J^ b {y). Substituting this into eq. (|3.2| ) and 
taking the limit p 2 — > we can rewrite r.h.s. of eq. ( p.2|) in the form: 

j- J d%(n a (P)\T{d,jt(y)4>(x) h T - ^(0)}|0) = 

/ d*y{8(y°)(Tr a (PMx) [J«\y), ^(0)]|0) (3.3) 

+6(y°-x°)(n a (P)\[J 5 b (y), j(x)) n y V(0)|0>}. 

The equal time commutators entering the above equation are simply the global axial 
transformations of the quark fields: 

Jd 4 y6(y )[J° b (y),m] = 4 75 ^ (0) ' 
f dS5{y°-x G )[4\y)^{x)] = #(*)y 7 5 • (3.4) 
Substituting this result into eq. ( |3.4| ) we get the following soft pion theorem: 

lim nOUt (-K a ( Pl )AP2) \$(x) V(0) |0> = ^(vr a (P)|^(x)n 75 ^(0)|0) . (3.5) 

Applying this soft pion theorem to the matrix elements entering the definition of the 
chirally even 27rDA eq. Q2.1|) and comparing the resulting expression with definition of the 
pion DA, eq. ( |3.1|) , we get the following relation between the pion DA and the isovector 
chirally even 27rDA: 

$f =1 (z,C = l,w 2 = 0) = -$f =1 (z,C = 0, W 2 = 0) = ip n (z) . (3.6) 
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The analogous theorem for the isoscalar part of the 27rDA's has the form: 

${=°( 0; ( = l, w 2 = 0) = ${=°(z, ( = 0, w 2 = 0) = . (3.7) 

Some comments are in order here. The derivation of the soft pion theorems presented 
here can be easily extended to other cases when we have an emission of soft pions. We 
should only keep in mind that for the case when we have nucleons in the out state (a 



typical example is a generalized skewed parton distributions, see [12, [I3|) the singularities 
related to the emission of soft pions from the nucleon legs must be taken into account. 
The soft pion theorem eq. ( |3.6|) allows us to relate the coefficients B^w 2 ) (see eq. fl2.9|) ) 



and the coefficients in the expansion of the pion DA in the Gegenbauer polynomials 

^(z) = Qz(l-z)(l + £ a^CT{2z-l)). (3.8) 

^ n=even ' 

The relation has the form: 



a 



W = £i?!S' =1) (0)- (3-9) 



l=i 



The soft pion theorem for the isoscalar chirally even 27rDA eq. fl3.7| ) implies that 



n+l 



EC u, (o) = o. (3.10) 

1=0 

4 Crossing relations between 27rDA's and quark dis- 
tributions in a pion 

An interesting feature of the pion is the fact that the 27rDA's are related to the skewed 
parton distribution in a pion by crossing. On one hand, the crossing relations derived in 
this section [] allow to use the information about quark distributions in a pion to constrain 
two-pion distribution amplitudes. On other hand, the measurements of 27rDA's in hard 
exclusive reactions can be used as a source of information on skewed quark distributions 
in a pion. Additionally the crossing relation are helpfull in modelling of skewed parton 
distributions. 

The matrix element entering in the definition of 27rDA fl2.1|) by crossing is related to 
the matrix element in definition of skewed distribution: 

^exp[iXrn ■ (p + p%7r a (p')\T{^ f (~Xn/2)h^ f (Xn/2)}\Ap)} = (4-1) 
5 a b 5 ff ,H I=0 (T, e, t) + ie Mc r c ff ,H I=1 (r, £, t) , 

where £ is a skewedness parameter defined as: £ = —(n-(p' — p))/(n-(p'+p)),t= (p'—p) 2 
and light-cone vector n M introduced in section |^. 

3 Part of the material of this section was presented already in recent preprints pi, Eq|. 
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This allows to express the moments of skewed parton distributions in terms of expan- 
sion coefficients B nt in eq. (|2.9| ): 

/ dTT N -i H \T,Z,t)=Y,Y,BUm N Cl /2 (l) / dx{\l-x*]x N ^Cr{x). (4.2) 

J - 1 n=0 1=0 VC/ ^-1 4 

Let us note that this expression satisfy the polynimiality condition for moments of SPD's 
(see ref. 0). 

In order to prove the relation (|4.2| ) let us consider the operator expression for the 
N— th moments of SPD. It is given by off-forward matrix element of local operator of spin 
N and twist-2: 



j^drr^H^t) = [( ^ + p )+] > W(V>-Vb) , (4.3) 

Analogously the iV— th moments of 27rDA ( |2.1| ) is expressed in terms of the following 
matrix element of the same local operator: 



dz(2z-l) N - 1 ^ 1 (z,C,w 



( P 'j#7 + (v r-vio> 



■ Y_J " + 1 .2^1/2 (W-PY 



n=0 1=0 



{p' +p) + 



\dx\[l-x"]x N -'cT{x), (4.4) 



where in the last equality we used the double expansion (|2.9|) . The matrix elements in 
eq. (|4.4j) and in eq. (|4.3|) are related to each other by crossing: 



W-p) + ] N ^ f B i {{v , _ v?)c i/i ( (p' + p) + 



n=0 (=0 

This is exactly relation ( f4.2|) . Let us note that in the expression ( |4.2|) the coefficients B n i(t) 
enter at negative argument whereas in ( |2.9| ) at positive argument. The corresponding 
analytical continuation can be done with help of dispersion relations (|5.8|) (see the next 
section). 

If we take the forward limit in (|4.2|) , we obtain the relations between moments of quark 
distributions in a pion and coefficients B n i(0): 

MP = f 1 dxx N - l (q 7r (x)-q n (x))=B I N =1 1N (0)A N for oddiV, 
Jo 

■ 2/ dx x N ~\q n {x) +qJx)) = Bir° 1N (0)A N for even N, (4.6) 
Jo 
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5/3, etc.) 
M 1 ( " ) = 1 

what corresponds to normalization condition (2.8). Also it is easy to see that Sf^°(0) = 
5/9M^ corresponds to normalization for isoscalar 27rDA []: 



where are numerical coefficients {e.g. A\ = l,A 2 = 9/5, A s = 6/7, A^ 
and q 7T (x) = m 71 " 4 (x). For example, from eq. ( |4.6| ) we obtain that i?Qf 1 (0) 



f 1 dz(2z - l)$ /=0 (z, C, w 2 = 0) = -2 AfJC(l - , (4-7) 

where MJ" is a momentum fraction carried by quarks in a pion. 

Using relations (|4.6|) and soft pion theorem (|3.6|) we can express the coefficient £?2i(0) 
[it describes the deviation of isovector 27rDA from its asymptotic form C,w 2 ) = 

6z(l - z)(2C - l)F^- m -(u; 2 )] as: 



jB ^(0)=a^ ) --A^ 



(tt) 



(4.t 



From the expression for the moments of SPD ( |4.2|) we can formally construct the SPD 
itself as series: 



tf 7 (T,e,*) = E E 7(1 - ^(1*1 - lei) ^(t)c^ 2 (^)c z 1/2 (-) , (4.9) 



n=0 «=0 



where each term in the sum has a support — £ < r < £. Let us note however that it 
does not imply that the SPD H(t,£, t) has the same support because generically the sum 
( f4.9|) is divergent at fixed r and £. In the formal solution (|4.9|) each individual term of 
the series gives the contribution to the amplitude which behaves as l/£ z at small £ what 
would imply the violation of unitarity for the hard exclusive reactions. This "disaster" is 
avoided through the fact the series is divergent and one can not calculate the asymptotic 
behaviour term by term. This situation is similar to duality idea: the entire skewed and 
usual quark distributions are given by infinite sum over t-channel resonances. Of course, 
there might be "non-dual" or "two component" way to compensate the bad behaviour of 
the resonances by some equally bad behaviour of non-resonant background. Such "two 
component" model can be written as: 

H'^ct) = E E i(i4MM-^i)« /2 (^; /2 (j) 

|=0n=I+l S S S 

+ -^non-resonant ( T 5 £ ) ^) ) (4-10) 

where Lq is the maximal spin of resonances which are taken into account and 
-^non-resonant( r 5 £5 t) * s a non-resonant background. The latter can be modelled in terms 



of double distributions ||16|| , however if Lq > the non-resonant contribution has to be 
oscillating function in order to avoid double counting in calculations of low moments of 
SPD (N < Lq) which are saturated by resonances in t-channel. 

4 To see this one has to use additionally soft pion theorem (p/fl) for isoscalar 27rDA 
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To summarize, the representation of SPD's ( f4.9| ) in terms of only resonance exchages 
gives a new picture of quark distributions in a pion which is dual to the picture based 



on the concept of double distributions [TO. It would be extremely interesting to build a 
model where the duality relation 

oo n+1 o 2 i 

,4 ^ U n=0 1=0 4 S S S 

is realized explicitly. 

At the end let us note that from the formal solution ( |4.9|) and soft pion theorems 
( |3.6| , |3.7|) we can determine the shape of SPD's in a pion in a limiting case: £ = 1 and 
t = 0: 

tf' =1 (r, £ = 1, f = 0) = i^(^yi) , # /= V, £ = M = 0) = , 

where ^^(z) is a pion distribution amplitude. These limiting relations impose strong 
constraints on the models of skewed parton distributions. For example it is easy to see 
that the model based on double distributions proposed in |L6| does not satisfy these 
constraints. 



5 Two-pion distribution amplitudes in the resonance 
region and DA's of resonances 

Generally speaking, the 2n DA's are complex functions owing to the strong interaction 
of the produced pions. Above the two-pion threshold w 2 = 4m 2 the 27rDA's develop an 
imaginary part corresponding to the contribution of on-shell intermediate states (2ir, 4it, 
etc.). In the region w 2 < 16m 2 the imaginary part is related to the pion-pion scattering 
amplitude by the Watson theorem [ 17 1 . This relation can be written in the following form 
(in this section we omit the indices || and _L but show the isospin index, since all formulae 
in this section apply to both chirally even and odd 27rDA's): 

lm^(z,C,w 2 ) = J d(phase spacers, w 2 )* A^iC MC, P±) , (5.1) 

where v4^. 7r (£ / ,p'j_|^,pj_) is pion-pion scattering amplitude with isospin / in the s-channel, 
(' and p' ± are the momenta of the intermediate pion, the phase space integral in these 
variables has the form: 



/ <i(phase space) = / d(' [ d9±, (5.2) 

J 6471"^ Jo Jo 

where 8± is the angle between p± and p' ± . In order to find the form of the irn scattering 
amplitude in the light-cone variables we start with the usual partial wave expansion of 
the pion-pion scattering amplitude: 
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4„ = 327r^(2/ + 1) sin[ ( 5/( W 2 )]e^^ 2 )p i [cos(^ m )] , (5.3) 
i 

where 8 cm is the scattering angle in the centre of mass system, it can be easily expressed 
in terms of light-cone variables. The relation has the form (here for a massless pion for 
simplicity )Q: 

cos(9 cm ) = 1 - 2(C + C - 2C0 - VC(1-C)C'(1-C0 cos(6 ± ) . (5.4) 

Substituting the partial wave expansion ( |5.3| ) into eq. ( |5.1| ) and computing the integral 
over 9± using the identity: 



d6 ± P t [l - 2(C + C' - 2CO - VC(1 - 0C(1 - CO cos(^ 

2tt C} /2 {2( - 1)C/ /2 (2C' - 1) , (5.5) 
we get the following expression for the imaginary part of the 27rDA: 

lm$ I (zX,w 2 +iO) = V(2/ + l) C dC^\zX\w 2 TM8l(w 2 )]e i5 l {w2) 

x C, 1/2 (2C-1)Q 1/2 (2C / -1). (5.6) 
If we now substitute the expansion of the 27rDA in the Gegenbauer polynomials (|2.9| ) in 



eq. (|5.6|) we get the equation for the imaginary part of the expansion coefficients B, 



ni ■ 



Im B^w 2 ) = sm[d( (w^je^B^w 2 )* = ^(w^ReB^w 2 ) . (5.7) 

One sees immediately that the phases of the coefficients B^ are the same as the phase 
shifts 5f . Let us also note that the equation ( |5.7|) for the case n — 0, 1 — 1 and I = 1 is 
identical to that obtained in |jTB[ for the pion electromagnetic form factor, which is not 
surprising because of identity B i(w 2 ) = F n (w 2 ). Since the form of equation (|5.7|) is very 
similar to that for e.m. pion form factor one can analyze it following ref. 

Using eq. ( f>.7\) one can write an iV-subtracted dispersion relation for the B^w 2 ): 



B^w 2 ) = £ 



N-l 2k (Ik 2N 

si(o) + " 



fe=0 



k\ dw 2k 



7T 



- ^ tan#(s)Re£?^(s) 

4m| sN ( s - W 2 - lO) 



(5.1 



The solution of such type of dispersion relation was found long ago by Omnes |T9[ and 
has the exponential form: 



^Expressions for massive case see in Appendix. 
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A> 2 ) = *<0)«p{g ^log *(0) ^ ^g-o) } ■ (5-9) 

This equation, by derivation, is valid below the inelastic threshold w 2 < 16m 2 , however 
the inelasticities (at least for the isospin one channel) in the pion-pion scattering are small 
up to rather large energies. Also the equation ( |5.9p being applied to the case of the pion 



form factor (n = 0, I = 1) |18|, 20] gives an excellent description of experimental data up 
to w 2 ~ 2.5 GeV 2 with two subtractions which are fixed by the value of the pion e.m. form 
factor at zero and its charge radius (see figure in p0|). Moreover the solution (|5.9|) can be 
systematically improved by taking into account higher thresholds (see discussion in |[21||), 
probably the most important is the KK one. The contribution of higher intermediate 
states can be suppressed by choosing sufficiently large number of subtractions N, see 
discussion in section |5TT| . 

A great advantage of the solution eq. ( |5.9|) is that it gives the 27rDA's in a wide range 
of energies in terms of known tttt phase shifts and a few subtraction constants (usually two 
is sufficient). The key observation is that these constants are related to the low-energy 
behaviour at w 2 — > of the 27rDA's. In the low energy region they can be computed 
using the effective chiral quark-pion lagrangian. In the next section we shall compute the 
subtraction constants using the low-energy effective lagrangian derived from the instanton 



model of QCD vacuum . For the first estimates of the 27rDA in this model in low energy 
region see ||. After the subtraction constants are fixed (see section ||) we shall know the 
27rDA's in a broad region of invariant masses. Of particular interest is the case when w 2 
is close to the mass of some resonance (say p): in this case we shall be able to extract DA 
of this resonance from the solution (|5.9|), this is explained in the next subsection. 



5.1 Distribution amplitudes of resonances 

The light-cone distribution amplitudes of mesons are defined as matrix elements of quark- 



antiquark non-local operators at light-like separations ||23|| . Here we consider as an ex- 
ample the p-mesons DA's, all equations in this section can be easily generalized to other 
mesons. Specifically, we consider chirally even and odd leading twist p-meson DA's defined 
as: 



(p\P)\j>{x)n T - V(0) |0) 



, + ™^=(»- (A, )^/>^ ,p '-v;w, (5.io) 



T 



<p°(P)|#r)^n^-#))|0> 



= i ( n . e W)i£L± dze iz ^^Az), (5.11) 



where is the polarization vector of the p-meson. The normalization constants and 
/J have dimension of mass and are chosen in such a way that (below we suppress indices 
_L and || where it does not lead to confusion): 
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dzcp p (z) = 1 . 

The light-cone p meson DA's are used to describe hard leptoproduction of vector 



o 



mesons 



24|, and in applications of light-cone QCD sum rules for exclusive semileptonic 



and radiative weak decays of the 5-meson p5[ 



There are certain problems with the definitions eqs. ( |5. lOlfTTTl) , because the p-meson 
is an unstable particle. One can treat the p-meson as a quasistable particle, i.e. consider 
the formal limit V p —>■ 0. From experimental point of view p- meson is detected as a Breit- 
Wigner peak in the invariant mass distribution of produced pions. Therefore a more 
adequate (and accurate) way to describe the hard production experiments is in the terms 
of two-pion DA's. 

In order to extract the rho meson DA's from the 27rDA's we use the expression (|5.9| ) 
of 27rDA in terms of the pion-pion phase shifts. In the vicinity of the resonance the 
corresponding phase shift (we shall speak for definiteness about p-meson, i.e. I = 1, 
isospin one) can be approximated by the Breit-Wigner formula: 



in r 

Sl(w 2 ~ ml) = arctan( / p ) . (5.12) 



m 2 — w 2 ' 



In the limit of vanishing width the contribution of the resonance part of the phase shift 
in the integral in eq. (|5.9|) can be written as : 



w 2N r°° , 6{{s) 



lim / ds— — = %- - log(m - w ) + logm (5.13) 

+ (pol. of order N — 1) . 

This result being exponentiated gives rise to the meson pole. Matching the solution for 
27rDA's eq. (|5.9|) near w 2 ~ m 2 to the Breit-Wigner form we get the a relation between 
p meson DA's and 27rDA's. We obtain the following expression for the coefficients of the 
expansion of p meson DA'S in Gegenbauer polynomials 



<p p (z)=6z(l-z)[l+ £ a^C 3 J\2z-r 

^ n=even 

in terms of subtraction constants entering eq. ( |5.9|) only 



N~l 

a(f) = 5 nl (0)exp(Eci nl) mf)' ( 5 - 14 ) 

k=l 

where the subtraction constants cj? 1 ^ can be expressed in terms of B n i(w 2 ) at low energies 

(5.15) 



w 2 =0 
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The constants /J' 1 - can be computed as: 

J o 



V2T p lmBl(m 2 p ) 



9pn 



# = V " P P ? 1V '^ , (5.16) 

9ptvttJ2tt 



where g p7rn is a p meson decay constant into two pions, is defined by eq. (£73] ). The 
constants B^{m 2 p ) are purely imaginary and can be computed using eq. ( |5.9| ) in terms of 
the pion-pion phase shift 8\ and low-energy subtraction constants. The latter are actually 
pion charge radius and radius of pion tensor form factor (they are estimated in the next 
section). 

eqs. ( |5.14| , |5TT6f ) can be easily generalized to the DA's of higher spin resonances. For 
example, the chirally even DA's of an isovector resonance of spin j can be computed in 
terms of the 27rDA's as: 



fc(z) = 6z(l - z){Cf_\{2z - 1) + E ^C 3 n /2 (2z - 1) ) , (5.17) 

v n=j+l 



with coefficients a$ computed as: 



(j) _ B nj (0) (Sr t (nj) 2k 



«■■■ — rexp(^4 n X)' (5.1J 
k=i 



The corresponding normalization constant is: 

f „ = ^^a=iiM) . (5 . 19 ) 

QRtvk 

The scale dependence of fn and 0$ is obvious. 

We see that the p meson DA's (we speak for definiteness about p meson, but all 
arguments below apply also to any other resonance) are fixed in terms of the subtrac- 
tion (low-energy) constants, which would be all fixed if we knew the exact low-enregy 
lagrangian of QCD0. The latter is not fully known, so we shall resort to the low-energy 



lagrangian derived from the instanton model of the QCD vacuum [22 



The r.h.s. of eq. Q5.14 ) depends on the number of subtractions iV whereas its l.h.s. is 



independent of this number. At first glance this fact implies that the subtraction constants 
are all zero. However we should keep in mind that the dispersion relation eq. (|5.8| ) 
is not exact, in its derivation the contribution of intermediate states other than two pion 
was neglected. In order to suppress the contribution of the higher states we have to take 



6 For example, the constant cf*' 01 ^ is proportional to the parameter Lg of Gasser-Leutwyler parametriza- 
tion of the EChL |30|. 
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a large enough number of subtractions in the dispersion relation (575) . This implies that 



one can expect an (approximate) independence of the r.h.s. of eq. (|5.14 ) of N only for 
sufficiently large N. 

Increasing the number of subtractions suppresses the high energy region, where the 
approximations were made, and hence increases the accuracy of the approximate solution 
( |5.9| ). However, increasing the number of subtractions we emphasize the low energy region 
the quantitative description of which is limited by incomplete knowledge of the effective 
chiral lagrangian. Therefore for phenomenological applications one has to choose some 
optimal N. The analysis of the pion e.m. form factor shows |L8], ^Of that already for 
N = 2 the contribution of the high energy states is sufficiently suppressed to describe the 
pion form factor with good accuracy up to w 2 ~ 2.5 GeV 2 . 

In the next section we take N = 2. Since with this value the contribution of the high 
energy states in dispersion relation (5.8) is sufficiently well damped and simultaneously 
we can reliably use the approximate low energy quark-pion chiral lagrangian derived from 
the instanton model of QCD vacuum ||22|| . In language of the chiral perturbation theory 
N = 2 corresponds to the fixing of the effective chiral lagrangian (EChL) to the fourth 
order. To this order the effective chiral theory derived from the instanton model gives 
22] , |26f the constants of the fourth order EChL (actually only its part describing pions 



interacting with external vector and axial currents) which are very close to ones obtained 
in the phenomenological analysis of Gasser and Leutwyler 



30 



6 Fixing low-energy constants from instantons 

The technique how to compute matrix elements of bilocal quark operators at low energies 
in the effective low-energy quark-pion theory was explained in details in refs. [27, ^ (for 
pion and photon DA's) and in || (for 27r DA's). Here we only sketch the main points of 
the calculations, for details the reader should consult refs. |27J, H. 

We shall compute the 27rDA's at small w 2 , eq. (|2.1| , |2.2|) , using the effective low-energy 
theory of pions interacting with massive "constituent" quarks, which has been derived 



from the instanton model of the QCD vacuum |2"2"|. The coupling of the pion field to the 



quarks, whose form is restricted by the chiral invariance, is described by the action 



S int = / d 4 x tfj(x)^M{d 2 ) W 5 (x) ^M{d 2 )^j{x), (6.1) 



where 

U^(x) = e ^r^(x)//. = x + * K a (x)T a - ir a (x)ir a (x) + .... (6.2) 

Here, f n = 93 MeV is the weak pion decay constant 

The momentum-dependent dynamical quark mass, M(p), plays the role of an UV 
regulator. Its form for Euclidean momenta was derived in ref . ||22|| . The momentum 
dependent mass cuts loop integrals at momenta of order of the inverse average instanton 
size, fT l ~ 600 MeV. One should note that the value of the mass at zero momentum, 
M(0), is parametrically small compared to p~ l ] the product M(0)p is proportional to the 
packing fraction of the instantons in the medium, which is a small algebraic parameter 
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fundamental to this picture. Numerically, a value M(0) = 350 MeV was obtained in 
ref. |2|. 



6.1 Isovector 27rDA's 

Using the technique described in M we can easily see that the loop integral in the effective 
theory for the chirally odd isovector 27rDA (|2.2| ) contains no divergencies, and hence we 
can neglect during its computation the momentum dependence of the constituent quark 
mass. Chirally even isovector 27rDA ( |2.1| ) was computed in || , the computation there was 
done neglecting the momentum dependence of the constituent quark mass and absorbing 
the divergencies into f n . It is easy to see that all (logarithmic) divergencies are contained 
in the w 2 -independent piece of this DA, the w 2 dependent piece is finite. The result of || 
in this approximation for w 2 = was very simple: 



$\\(z,(,w 2 = 0) = sign(l - 2C)0[(z - C)(l - C - z)\ , (6.3) 

where 6(x) is a step function. We see that chirally even isovector 27rDA exhibits jumps at 
points z = C and z — 1 — £. As was shown in |J7|, |27|, |2S] the momentum dependence of the 
constituent quark mass becomes crucial at these points. Here for the numerical estimates 
of the w 2 independent piece we employ a simple numerical fit to the momentum dependent 
constituent quark mass obtained from the instanton model of the QCD vacuum, 

(1 + 0.5p 2 p z ) 

As we already mentioned the ^-dependent piece of the chirally even isovector 27rDA 
is finite (it does not exhibit the jumps) and hence it can be computed neglecting the 
UV-cutoff provided by the momentum dependence of the constituent quark mass. 

Results of a computation of the chirally even isovector 27rDA at w 2 = are shown in 
Fig. [TJ, where we have plotted $| =1 (z, w 2 = 0) as a function of z at various values of 
(. We see that the function is concentrated in the interval of z between ( and 1 — ( and 
the jumps are smoothed. Also it is interesting to note that if we take ( = 1 the result for 
$| =1 (z, (, w 2 = 0) coincides exactly with that for the pion DA obtained in |27, 28] in the 
same model, what is not surprising due to the soft pion theorem ( p.6|) . 



Inspecting loop integrals for the chirally even 27rDA's in the model (see typical expres- 
sions in ref. |§) we can make interesting observation: in the limits £ — > 1 and — * the 
w 2 dependence appear exclusively in the combination z(l — z) w 2 . This observation allows 
us to obtain the following model relations for the coefficients of Gegenbauer expansion 
b\i{w 2 ) (we omit the superscript || below to make the expression not so busy): 

Y.^p j B nT 1 {w 2 )\ w 2 =0 = 0, for any even n > 2j , (6.5) 

odd 

]T— l^fV 2 )!^ = 0, for any odd n > 2j - 1 . (6.6) 



even 
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These model results, in a sense, can be viewed as a generalization of the soft pion theorems 
( |3.9| , |3.10|) . Although we find these relations in the model calculations, we can speculate 
that they have more general nature and follow from the chiral Ward identities in the limit 
of large colour number N c — > oo. This conjecture will be worked out in details elsewhere 

Below we give the results of numerical calculations of a few first coefficients of Gegen- 
bauer expansion ( ^.9| ) B^^w 2 ) of the isovector chirally even 27rDA at low w 2 : 



B ' kw2) = °' 26+ Tiift + "" (6J) 

II , 9x llN c w 2 

B! l(w >) = -0.006 - 5S5 ^ + ..., (6.10) 

1 1 N vi 2 

- -^-m£k + -- < 6 - u > 

ii / 9x HN c w 2 

Bkw) = °- 14 + 7560% + --- (6 - 12) 

From this calculation we immediately can extract the coefficients of the Gegenbauer ex- 
pansion of the pion DA (p.8|) using the soft pion theorem (K 



(tt) 
«2 



0.06, 



(tt) 
04 



0.014 . 



(6.13) 



These are exactly the values obtained in [28 . 

Remembering that B^w 2 ) = F 7T (w 2 ) we extract from eq. (6.7) the value of the pion 
electromagnetic charge radius (without chiral loops!) 



(r 2 ) e . m . 



(6.14) 



which coincides with the result obtained previously in ref. P2"| . The model result ( 6.14 ) 
numerically gives (r 2 ) e m . = 0.35 fm 2 to be compared with the experimental value 0.439 ± 
0.008 fm 2 ||29| . Comparing these numbers we should keep in mind that the model calcula- 
tions do not include the chiral logarithms, the latter appear when we perform the integral 
in eq. ( |5.9| ) (for n — 0, I — 1) with the pion-pion phase shift obtained in the leading order 
of the chiral perturbation theory: 



w 2 {l-Am 2 Jw 2 fl 2 
96nf 2 



Then the result for the pion charge radius becomes (remember that Bqi(w 2 ) = F- K {w 2 )): 
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1 



log( 



0.45 fm 2 



mi 



47T2/2 l 67r2/ 2 

in good agreement with the experimental value. 

Basing on the results ( |6.7| - |6.12|) we can easily estimate the low-energy subtraction 
constants, see eqs. (pTS|j5lgD , the results are summarized in Table 1, they refer to the 
normalization point p ~ l/p w 600 MeV the inverse instanton size in the instanton 
vacuum. 

Having the values of the low-energy subtraction constants we can obtain the chirally 
even isovector 27rDA in a wide interval of invariant two-pion masses using the solution 
( |5.9[ ) with N = 2 (see the discussion in previous section). The first qualitative observation 
we can make on basis of Table 1 is that the coefficients in front of C^ 2 (2z — 1) for n = 2, 4 
(they describe the deviation of DA's from its asymptotic form) for the P-wave production 
of pions are negative and B^w 2 ) increases (by absolute value) whereas B^iw 2 ) decreases 
with increasing of w 2 . 

Here we quote the result for the chirally even p meson DA extracted from 27rDA using 
eqs. ( CTJ5TT6D : 



<pUz) = 6z(l -z)[l- 0.14C 2 3/2 (2^ - 1) 



0.01 cT 



(2* -i; 



(6.15) 



with normalization constant /J = 190 MeV according to eq. ( [5 .16 ) in a good agreement 



with experimental value /j = 195 ± 7 MeV |3*I||. The p meson DA's were the subject 



of the QCD sum rules calculations |23|, |32|, |33], our result eq. ( 6.15 ) is in a qualitative 



disagreement with the results of QCD sum rule calculations, the sign of a\ obtained here 
is opposite to the QCD sum rules predictions a| = 0.18 ± 0.1 |32| and a| = 0.08 ± 0.02, 
a\ = —0.08 ± 0.03 |33] (these results refer to normalization point p = 1 GeV). The 



difference in the sign of ours and QCD sum rule calculations can be traced back to the 
signicicantly different values of . Indeed, using the eq. (|5.14| ) and the crossing relation 

and the third moment of the quark distribution in a pion in 



( |4.8|) we can relate a 2 
the following way: 



J B^ =1 (0)exp(cf 1) m 



i?r(0))exp(ci 21 W) 



(a^ - -M^) exp(cr m 2 p ) ■ 



(6.16) 



The value of m the instanton model of the QCD vacuum is very small (see eq. ( 6.13 ) 



and ref. p8[), such that c4 < 1-^3 ■ O n contrary the value of in the QCD sum rule 
calculations 0, |35| larger than ^M^\ what implies the positive sign of of . 

Using eq. ( |5.14| ) and the value of the subtraction constants from Table 1 we can make 
a rough estimate of the shape of the p' chirally even DA (m p i = 1.465 GeV): 



6*(l-z)[l 



0.055 Cl /2 {2z 



1)- 0.064 C 3 4 /2 {2z-l) + ...} 



(6.17) 
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This is just a rough estimate because the mass of p' is rather big and importance of 
neglected in the derivation of eq. ( |5.9|) high energy states increases and hence the accuracy 
of the solution eq. ( |5.9| ) becomes poor. But even from this rough estimate we see an 
interesting tendency: the coefficients in front of the Gegenbauer polynomials with n > 2 in 
the expansion of p' meson chirally even DA become comparable with the leading nontrivial 
(n = 2) coefficient and the DA develops nodes. 

Now we turn to results for the chirally odd 27rDA eq. ( |2.2j ). Simple calculations of 
matrix element eq. fl2.3| ) give the following value of the constant j-^ 

8n 2 f 2 

= -j^ « 650 MeV . (6.18) 
Calculations of B^(w 2 ) in the low energy region are straightforward (details will be 



published elsewhere PSfl), the results are 



B 01 


(w 2 ) 


itr 

+ 12M 2 + ' ' ' ' 


B 21 


(w 2 ) 


7 w 2 
36 ^ 30M 2 + ' ' 


B 23 


(w 2 ) 


36 ^ + 30M 2 + ' ' 


B ii 


(w 2 ) 


11 5w 2 
225^ 168M 2 + 


B 43 


(w 2 ) 


77 w 2 
675 ^ 630M 2 + 


B h> 


(w 2 ) 


11 w 2 
135^ + 56M 2 + ' 



(6.19) 
(6.20) 
(6.21) 

), (6.22) 
), (6.23) 

(6.24) 

From these equations we get immediately the values of the subtraction constants c^ nl \ 
they are summarized in Table 1. The value of "tensor" isovector radius of pion is from 
eq. (|6.19|) (without chiral loops!): 



(r 2 )t = ^p^ 0-16 fm 2 . (6.25) 

The chiral loops are the same as for charge radius. 

Using the results from Table 1 for c^ nl ^ and eqs. ( 5.14 , 5~16|) we obtain the chirally 
odd p-meson DA: 



<p£{z) =6z(l-z)[l + 0.11 Cl /2 (2z-l) + 0.03 C 4 3/2 (2^ - 1) + . . .] , (6.26) 
and the normalization constant fr: 



771 

fp =f l T £ e M((r 2 )t - (r 2 ) e . m .)m 2 /6] « 140 MeV . (6.27) 

J2vr 
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[This value refers to the normalization point /i ~ 1/p ~ 600 MeV] . The QCD sum rule 



calculations ||, |3[ give = 160 ± 10 MeV af = 0.2 ± 0.1 || and = 169 ± 5 MeV 
= 0.36 ± 0.03 [33| at the normalization point /i = 1 GeV in fair agreement with our 
results. 

The results for parameters of the p meson DA's are summarized in Table 2 where we 



also quote the results of the QCD sum rule calculations [32, p3| . With the help of Table 1 



and eqs. ( |5.9| , |5.18[ ) the reader can easily construct DA's of isovector mesons of any spin. 
6.2 Note on isoscalar 27rDA's 

The quark isoscalar chirally even 27rDA's under ERBL evolution || mix with gluon 
27rDA's. In the instanton vacuum model considered here the gluon 27rDA's are para- 
metrically small at low normalization point (p ~ 1/p) in the instanton packing fraction 
(see ref. f38j) and hence are zero at the level of accuracy considered here. 



The second moment of the chirally even 27rDA eq. ( |2.1|) is related to the quark part 
of the (symmetric) energy-momentum tensor T£ v : 

J o 1 dz(2z-l)^f°(z } Cw 2 ) = ^ 2 out(K a ( Pl )APi)\T q ++ \0) ■ (6.28) 

In the instanton model of the QCD vacuum (see |38[) the gluon part of the energy- 
momentum tensor (projected on the light-cone direction n M !) is parametrically small in 
the packing fraction of the instanton liquid. Therefore here we can use, the so-called, 
quark- ant iquark approximation, i.e. T+ + = T ++ + small corrections. At this level of 
accuracy we have (cf. equation (fl.7|)): 

£ dz(2 z-l) $j =0 (^, C, w 2 = 0) = -2C(1 - C) • (6-29) 

The model calculations of the isoscalar chirally even 27rDA give the following result 
for its second moment: 

rl NT 7/ i2 

/ dz(2 z-l)$J=°( 2 ,C, W 2 ) = -2C(l-0[l + T^7i + •••]• ( 6 - 3 °) 
Jo 11 487i - z j* 

We see that the overall normalization is fixed by the energy momentum conservation and 
the radius of pion form factor of the energy momentum tensor is: 

(Oemt = = 0-18 fm 2 . (6.31) 

From eq. ( |6.3U| ) we get immediately the values of the low-energy subtraction constant: 

4,(0) = — M' 2 (0) = | , (6.32) 
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A 1 °g S fo(«' 2 )L=»=o= AlogM' 2 (^ 2 )L 2= o= T^L = 0.73 GeV" 2 . (6.33) 
aw z aw z 48tt z j z 

Obviously the soft pion theorem eq. ( p.lO|) is satisfied. Let us use the above model 
calculations and crossing relation ( |4.2j ) to estimate the low momentum transfer behaviour 
of the second moments of skewed quark distribution in a pion. From eq. (|4.2|) for the 
second moment we have 

£ dxxH^x, e, t) = ^j (i) e + ^BUt) e c 2 l/2 Q 

= jj[4^H 2 + 4^)^]- (6-34) 

Now if we substitute in the above equation the result of the model calculation ( |6.33| ) we 
get: 



r 1 1 N t 

jjXXH^{X^t) = -(1 - e 2 )[l + ^ + ...]■ (6.35) 

7 Two-pion distribution amplitudes in hard pions pro- 
duction processes 

The dependence of the two-pion hard production amplitude ( |1.3|) on the 27rDA's factorizes 
into the factor: 

^oc /'-^Sf^C^O 2 ), (7.1) 

JO Z{L — z) 11 

for the two pions in the isovector state, and 



A oc /* dz^—^r Sf=°(z, C, w 2 ; Q 2 ) , (7.2) 
Jo z(l — z) 11 

for the pions in the isoscalar state (e.g. for the 7r°7r° production). In the above equations 
we showed also the dependence of the 27rDA's on the scale of the process Q 2 , which is 
governed by the ERBL evolution equation [fj. 

For the process ( p..3|) at small xsj (see e.g. recent measurements [p9| ) the production of 
two pions in the isoscalar channel is strongly suppressed relative to the isovector channel 
because the former is mediated by C-parity odd exchange. At asymptotically large Q 2 
one can use the asymptotic form of isovector 27rDA: 



lim <S> I f 1 (z,(,w 2 ;Q 2 ) = 6F n (w 2 )z(l-z)(2(-l) , (7.3) 

where F n (w 2 ) is the pion e.m. form factor in the time- like region. Therefore the shape 
of 7r + 7r~ mass spectrum in the hard electroproduction process ( |1.3| ) at small xbj and 
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asymptotically large Q 2 should be determined completely by the pion e.m. form factor in 
time-like region: 



lim A oc e iS ^ w2) \F n (w 2 )\(2C - 1) . (7.4) 

Q 2 ^oo 

Deviation of the 7r + 7r~ mass spectrum from its asymptotic form eq. (|7.4j) ("skewing") can 
be parametrized at small Xsj and large Q 2 in the form: 

A oc e^\F n (w 2 )\ l + fi 21 (0; /*) exp{ci 21 V} ( J ( 2 C - 1) + 

Va s (/io)/ J 

e^B 23 (0;p ) exp{^ 2 + ^(^)}(^T) 50/(99 " 6iV/) C3 1/2 (2C - 1) 

+higher powers of 1/ log(Q 2 ) , (7.5) 



here 



R\{w 2 ) = Re— / ds— — — — and 6 ni = log 



.2> 



/ •"■ o, > , ., ->-',-, "/;(■-" > • (7-6) 

We see that the deviation of the nir invariant mass spectrum from its asymptotic form 
eq. ( |7.4j ) in this approximation is characterized by a few low-energy constants (I?2i(0), 

(21) 

-^23(0), c\ , 623), other quantities- the pion e.m. form factor and the 7777 phase shifts-are 
known from low-energy experiments. In principle, using the parametrization (|7.6[) one 



can extract the values of these low-energy parameters from the shape of rnr spectrum 
(not absolute cross section!) in diffractive production experiments. Knowing them one 
can obtain the deviation of the n meson DA (see eq. (|3.9| )) 

= B 21 (0) + B 23 (0) , 

and the p meson DA (see eq. Q5.14D ) 

a 2 p) =B 21 (0)exp(cf ) m 2 p ), 

from their asymptotic form 6z(l — z), as well as the normalization constants for the DA 
of isovector resonances of spin three (see eq. ( |5.19| )). Additionally the parameter -823(0) 
is related by the crossing relations to the third moment of the valence quark distribution 
in a pion (see eq. (fO|)): 

M* = ^B 23 (0). (7.7) 

We see that, in principle, the diffractive production of two pions can be used to get an 
information about quark distribution in a pion. 

In analysis of experiments on two pion diffractive production off nucleon (see e.g. |39| ) 
the non-resonant background is described by Soding parametrization ET|, which takes into 



account rescattering of produced pions on final nucleon. Let us note however that in a 
case of hard (Q 2 — > 00) diffractive production the final state interaction of pions with 
residual nucleon is suppressed by powers of l/Q 2 relative to the leading twist amplitude. 
Here we proposed alternative leading-twist parametrization ( f7.6|) describing the so-called 
"skewing" of two pion spectrum. 
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8 Conclusions 



We showed that the two-pion distribution amplitudes are the most general object to 
describe the hard electroproduction of two pions. Using Watson final state interaction 
theorem and the soft pion theorems proved here, we can determine the 27rDA's in a wide 
region of pion invariant masses in term of the pion phase shifts and a few low-energy 
(subtraction) constants. The former are know from data on tttt scattering (soft physics), 
the latter are non-perturbative input which we estimated here using the instanton model 
of the QCD vacuum. We showed that these non-perturbative characteristic can, in prin- 
ciple, be extracted from the shape of nir mass spectra in diffract ive pions production 
experiments. Then they can be used to determine, say, the deviation of the tt meson DA 
(see eq. (|3.9|)) and p meson DA (see eq. ( |5.14|) ) from their asymptotic form 6z(l — z), 
and hence to obtain non-perturbative information about structure of mesons. We demon- 
strated (see section |5.1|) that DA's of resonances of any spin can be determined from 27r 
DA's. 

We derive the crossing relation which relate various distribution amplitudes to the 
quark distributions in a pion. 

The soft pion theorems considered in section |^ give an example of how low-energy 
(chiral) physics can be studied in hard processes (see recent work in this direction |40|). 

The methods developed here can be easily generalized for more complicated cases. One 
of examples is generalized N — ► ttN skewed parton distributions Hl2| , |13f entering the QCD 
description of the "non-diagonal" deeply-virtual Compton scattering + p — > 7 + ttN. 
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A Appendix 



Due to the kinematical constraint ( |2.5| ) the physical region of variable ( is bounded to 

the interval 

C- < C < C+, where 



1±\1- 



4m| \ 



w 



(a.i; 



It implies that the Gegenbauer polynomials Cj^ 2 (2( — 1) used in the expansion ( [2.9|) are 
not orthogonal on the interval £_ < ( < ( + . The correct set of polynomials in this case 
is: 



(A.2) 



It is easy to check that these polynomials are orthogonal on the interval < ( < (+: 



f + d( Vi{c)v k (0 



i^k. 



1 1 _ ±z4 7 _ I, 

22+1 V 1 7„2 5 



(A.3) 



In the case of w close to the threshold 4m^. the expansion ( |2.9|) is modified: 



oo n+1 



C, w 2 ) = 6z(l - z) £ E S^V)^** - 1)^(0- 

ra=0 Z=0 



(A.4) 



From the normalization condition 



we conclude that 



01 



^ m \pc.ui./„„2\ 



W 



2 7T 



[W 



(A.5) 



Also the equations in section [5] are modified near threshold. The phase space volume 



is: 



/ <i(phase space) = / rfC / ^J- • 



(A.6) 



The scattring angle in the cm. system is expressed in terms light cone variables as 
following 



cos 



1 — Am^/w 2 



i-2(c + C'-2CO 



- 4 V /(C(1 - C) - 2 )(C'(1 - CO - ™ 2 Jw 2 ) cos(^) 



(A.7) 
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With this equality the expression ( |5.5| ) is modified as follows: 



r d^[cos(# cm )] = 2ttP z (C) Vi(C) . (A.* 
Jo 



This demostrates that the set of polynomials ( |A.2j ) corresponds to the expansion in partial 
waves of produced pions. 
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(nZ) 


Bl(0) 


B^{0) 


c \\(ni) Ge y-2 


cf (nl} GeV" 2 


-A log Bl(0) GeV" 2 

a.w* <-> m \ / 


3^ log B^(0) GeV" 2 

a.w* a '« \ / 


(01) 


1 


1 








1.46 


0.68 


(21) 


-0.2 


0.2 


-0.6 


-0.95 


0.85 


-0.27 


(23) 


0.26 


0.2 








0.65 


0.27 


(41) 


-0.006 


0.05 


1.1 


-0.92 


2.5 


-0.24 


(43) 


-0.12 


0.11 


-0.36 


-0.28 


0.30 


-0.01 


(45) 


0.14 


0.08 








0.36 


0.15 



Table 1: Values of the subtraction constants entering eqs. ( ^9| , |5.14|) calculated in the 
model of instanton vacuum. Isovector case. 





our results 


Ball & Braun 


Bakulev & Mikhailov 


f\\ 

J p 


(MeV) 


190 


198 ±7 


201 ±5 


f ± 


(MeV) 


140 


160 ±10 


169 ±5 




al 


-0.14 


0.18 ±0.10 


0.08 ±0.02 




a\ 


-0.01 




-0.08 ±0.03 






0.11 


0.2 ±0.10 


0.36 ±0.03 






0.03 







Table 2: Parameters of the twist-2 vector meson distribution amplitudes extracted from 
the 27rDA's compared with QCD sum rule calculations [34,35]. Our results refer to nor- 
malization point fi « 600 MeV, QCD sum rule calculations to \i ~ 1 GeV. 
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2.8 




_Q 2 I ' I ' I ' I ' I ' 

' 0.0 0.2 0.4 0.6 0.8 1.0 



Figure 1: The two-pion chirally even isovector distribution amplitude <3>| 1 (z,(^,w 2 = 0) 
at various values of (. Solid line: ( = 0.1. Dashed line: ( = 0.3. Dotted line: ( = 0.45 
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